Unconventional Fulde-Ferrel-Larkin-Ovchinnikov states in spin-orbit coupled 

condensates: exact results 
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We find that a model Hamiltonian of s-wave superconductors in the presence of spin-orbit inter- 
actions and a Zeeman field is exactly solvable. Most intriguingly, based on the exact solutions, an 
unconventional type of Fulde-Ferrel-Larkin-Ovchinnikov (FFLO) ground state is rigorously revealed, 
in which the center-of-mass momentum of the fermion pair is proportional to the Zeeman field. We 
also elaborate on the drifting eff'ect of the Zeeman field on the spin-orbit-coupled Bose-Einstein 
condensate. 
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Topological superconductors arc expected to be po- 
tential platforms for topological quantum computation 
and therefore have attracted great interests recently in 
the fields of condensed matter, cold atoms, and quan- 
tum computation. As is known, gapless edge-modes are 
topologically protected against any weak disorder that 
cannot destroy the bulk gap and break the relevant sym- 
metries. It was long ago proposed that non-Abclian Ma- 
jorana fermions supported by the topological supercon- 
ducting phase may be realized in the chiral px + ipy su- 
perconductors [ij . Recently, a rather simple realization of 
topological superconductors was theoretically suggested 
by inducing the effective chiral p-wave superconductivity 
on the surfaces of topological insulators Q or semicon- 
ductor films Q in proximity to an s-wave superconduc- 
tor. Both suggestions are in the framework of mean-field 
theory (MFT) and involve crucially spin-orbit (SO) in- 
teractions related to Dirac fermions. Very recently, the 
pairing Hamiltonian of Dirac fermions in two dimension 
(2D) was shown to be exactly solvable Q and an equiva- 
lent chiral p-wave ground state was rigorously confirmed 
beyond the mean-field approximation. 

On the other hand, 3D topological superconductors 
[^-0] were very recently reported in the copper-doped 
topological insulator Bi2Se3. The Fermi energy of the 
superconducting sample was found to be in the relativis- 
tic regime of bulk band dispersion, indicating that the 
electrons participating superconductivity are also Dirac 
fermions [6|. Moreover, it was revealed that the effective 
spin-orbit interactions or Dirac fermions can also be sim- 
ulated with controllable experimental parameters using 
ultracold atoms in optical lattices both theoretically 
and experimentally [lo| , including 3D cases even with a 



rather flat band ll| . Therefore, it is natural and critical 



to ask: (i) whether the BCS-type model of Dirac fermions 
is exactly solvable for 3D cases; (ii) how to solve it if the 
answer is yes; (iii) whether the existing important con- 
clusions for the 2D are valid or not for the 3D; and (iv) 
whether an exotic quantum state of Dirac fermions can 



be emergent in the presence of Zeeman field. 

In the present work, we answer the above inter- 
esting and important questions timely and unambigu- 
ously. We study the SO-coupled superfluidity in the 
presence of an effective Zeeman fleld. Conventional s- 
wave superconductors subject to a sufficient strong Zee- 
man field might exhibit the so-called Fulde-Ferrel-Larkin- 
Ovchinnikov (FFLO) state [l2j. However, the conven- 
tional FFLO state is rather weak, such that concrete ex- 
perimental verifications of this intriguing state are still 
highly awaited. While, the situation is distinctly differ- 
ent for pairing Dirac fermions because the Zeeman field 
acts on the Dirac fermions just like an effective vector 
potential 0, , and consequently a new type of uncon- 
ventional FFLO state may be emergent, with the center- 
of-mass momentum of Cooper pairs being proportional 
to the Zeeman field. Here we solve the FFLO ground 
state exactly, and elaborate that it is stable at least for 
a special case and may be observed more easily in exper- 
iments. 

We consider a fermionic pairing Hamiltonian with spin- 
orbit interactions in a 3D system in the presence of a 
general effective Zeeman field B = {B^, By, B^) coupled 
with the spin degrees of freedom. The total Hamiltonian 
reads 



H — Ho + Hint , 



with 



^" = ^(4t'4i)(£k + [ak + AB] • a)(ckt, c^if, 

k 

-ffint = ^0(k,k')4+q^cLk+q4,C-k'+q4Ck'+qt, 
k,k',q 



(1) 



(2) 



where k = {kx,ky,kz) is the wave- vector and cr = 
{ax,(yy,Uz) with (Jx,y,z the pauli matrices. Ek = £-k 
is the spin-independent single-particle energy, a and A 
parametrize the SO and Zeeman coupling strengthes, re- 
spectively. Vb(k, k') = Vb < 0. For physical relevance. 



2 



the Zeeman field B may correspond to the laser-induced 
effective Zeeman field (9|, [14 1 in ultracold fermionic sys- 



tems (or the magnetic field applied within the plane of 
2D electron pairing systems). The single-particle part of 
the Hamiltonian, i.e. Hq, can be diagonalized by the uni- 
tary transformation (ckfjCkj.)"^ = J7k-i-b(ak+, Ok-)'^ with 
b = Aa~^B and 



cos t/k 
sin6'ke*'^ 



sm f^ke 
cos 0k 



(3) 



Here ^k S [0,7r/2], y^k G [0,27r). The definitions of 61k 
and (pk are given by 



cos(26'i 




k = 

otherwise ' 



and 



1, k. 







{kx + iky)/\kx -f iky\, otherwise 



(4) 



(5) 



After the above transformation of local spin rotation in 



the k-space, i/o = Z]k.s=± ^k,s 



£k s^-kiS is diagonalized with 



'k,s 



ek + salk + b| 



(6) 



From the above equation, it is clearly seen that for con- 



stant or weak k-dependent £k £ 



b 

k-b.s 



-k-b.s 



, SO that 



the intra-branch pairing (k — b, s; — k — b, s) with the to- 
tal momentum —2b is favored in the presence of the ef- 
fective Zeeman field and the spin-orbit interaction. The 
idea is schematically illustrated in Fig.[T] For general ek, 
the nesting vector q, satisfying Ek+q s ~ ^-k-i-qs' "^^^ 
equal to — b but still proportional to the Zeeman field 
Thus with the nesting vector q. the reduced BCS 




FIG. 1. Illustration of the evolution of the Fermi surface 
and energy shell with B for the case where Ek is constant 
or weakly k dependent. The shaded areas denote the energy 
shells around the Fermi surfaces. 

Hamiltonian is singled out from a gro up of candidates, 
which is related to the FFLO state [ll [H [H, 



-ffint(q) = V'o^C^qt'^-k+qiC-k'+qiCk'+qt. 
k,k' 



(7) 



The above equation can be rewritten in terms of ^ and 
flk,s by applying the unitary transformation and accord- 
ingly the pairing Hamiltonian becomes 

Hp = Ho + Vo J2 e-"'"''+'^'^'''4,,(qMk',s'(q), (8) 

ks,k's' 



where Hq is re-expressed as 

^0= E 



'k,s 



b,q - q 



k,s 



(9) 



k,s=0,± 



with 



-k+q 



k-|-q,s 



^Wk^ = { £^+q,+ + (-)£^+q,- 



and 



^k+q,s + ( — )^-k+q,s, 
%+q,+ + (-)"-k+q,-, 

0, 



k G A'l, s = ± 

k e A'o,s = 
otherwise 

k e i4:i,s = ± 
k e ii:o,s = . 

otherwise 



In the above formulas, i^o = {^\{kxiky^kz) ~ 0} and 



A'l = {k|fc^ 



0,fc^ > 0} U {k|fc^ = 0,fc.y > 0} U 



{kjfcj; > 0}. The pair operator in Eq. ([8]) is defined as, 
I+q,,aLk+q,., k e Ki and s = ± 
4,.(q) = { aUq,+«-k+q,-. k e /vo and .s = • 



0, 



otherwise 



(10) 

In the derivation of Eq. ([8]), we have used the relations 
cos6l_k = sin6'k and e"^-'' = -e"^'' for (fc^,fcy) 7^ 0. 
Note that the definitions of the above pair operators are 
slightly different from those in the 2D case which 
plays a crucial role in solving the 3D model exactly. 

The operators involved in the pairing Hamiltonian sat- 
isfy the following commutation relations 

Ks,4'..'(q)] = 25k..k'.'4,.(q)> 

[Ak,.(q),4',.'(q)] = W.s'(i-^L)' (11) 
["^k.^'^^k',.'] = ["^k,s'^k%s'(q)] = 0, 

and therefore the commutator algebra is closed. Using a 
scenario similar to that employed in Ref. the eigen- 
states of the Hamiltonian Eq. ([5]) are found to be 



|n,5+,5'_ 



n n 4,+q,-n^iio) 

(12) 



with 



B, 



E 

,s=0.±; 

keP^ 



2Ck\" 



(13) 
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where the parameters E^, 



1,2,.. 



are solutions 



of the ri-coupled Richardson's equations [17 



s=0,±: ^"^k.s 
kG-Pa 



Vo 



E 



2T^o 



En, — E^ 



(14) 



Here S± denotes the set of singly occupied levels (namely 
blocked levels) of the ± branch with cardinality to±, 
while P± and Pq the sets of levels with the blocked ones 
excluded. The state vector defined in Eq. (fT2|l describes 
the eigenstates of TVg = m+ + ?7i_ + 2n fermions with n 
as the number of pairs. The corresponding energy is 

n 

E{n,S+,S.)^ ^ ^k+q,+ + E ^k+q,-+E^- 



keS4 



keS- 



(15) 

For the special case of B = 0, it is straightforward to 
conclude that the SO-coupled 3D pairing fermion Hamil- 
tonian is exactly solvable, as in the 2D case Several 
remarks on the present 3D case can be made, (i) We 
calculate the gap function (as well as the pair-pair corre- 
lation function) as we did in the 2D case, and find that 
Ak,s = e"**'^'' where ^ is always real and just the 
usual s-wave one. Therefore the pairing symmetry of the 
gap function Ak,.s is of px + ipy-wa,ve. (ii) We can indeed 
show from Eq. (|12|) that the pairing ground state has the 
time-reversal symmetry, as expected by the MFT. (iii) 
The results can be readily generalized to the cases with 
anisotropic spin-orbit interactions {axk^, ayky , azkz) ■ cr 
by a scale transformation of k. (iv) If — 0, we are 
able to recover those results obtained in our previous 
2D work (v) It was indicated by Gaudin that the 
Richardson's equations ([M]) can recover the BCS gap 
and number equations in the continuum limit under the 
single-arc assumption [3, [3l • 

Most interestingly, when Sk is a k-independent con- 
stant, the Hamiltonian for B 7^ can be related to that 
for B = by the following simple mappings: 



flk-b,; 



A. 



4..(b) 



4,+rflk,- 
0, 



k e Ki and s = ± 
k e A'o and s = 
otherwise 



for the single-particle and pair operators, respectively. 
From this exact mapping, we have for the ground state 
with n pairs and no unpaired fermions. 



E,{B) = E,(B = 0), v = l,2,- 
|n)FF(B)=T_b|n)(B = 0), 
E{n,B) = E{n,B = 0), 



(16) 



where T_b means the translation of each fermion pair 
by —2b. The peculiarity of this type of FF state is 



that the center-of-mass momentum of the Cooper pair 
is Q = — 2b with its order parameter being proportional 
to e**^'^. This drifting effect of the Zeeman field on the 
superfluid of SO-coupled (or Dirac-type) fermions is sim- 
ilar to that of the vector potential of the magnetic field 
on a superconductor, but without having to increase its 
kinetic energy. And due to the nature of Dirac fermions, 
the mechanism of the above FF state is definitely differ- 
ent from that accounting for conventional FF states. For 
the isotropic SO-interaction, Q align in the same direc- 
tion of the Zeeman field, while for the anisotropic case, 
Q is along the direction of the scaled Zeeman field. 

As a simple example, we consider a special case where 
all N fermions are on the Fermi surface |k-|-b| = kp, 
which may be viewed as an approximation for the consid- 
ered system as many physical phenomena may be mainly 
related to the electrons near the Fermi surface. The de- 
generacy of the Fermi level is assumed to be ft. In the 
presence of the SO-interaction, the energy level is split 
into two branches ± ~ £_f ± akp. The symbols in 
Fig. [2] denote the scaled condensation energy AE/fl as a 
function of the fermion density for two values of and 
the scaled interaction fixed, calculated from Eqs. 
and numerically. These exact results of AE/^l are 
always higher than those obtained from the MFT and ap- 
proach to them in the large limit, as shown in the fig- 
ure. Moreover, some results are beyond the expectation 
of the MFT. We find that the exact condensation energy 
AE is always finite at half filling {6 = 1), indicating a su- 
perconducting state in contrast to the normal-state (for 
VqQ < 2) predicted by the MFT. The exact data is obvi- 
ously asymmetric with respect to (5 = 1, while the MFT 
curve displays a clear particle-hole symmetry. These dif- 
ferences between the exact solution and the MFT are 
more notable for small N cases and therefore more sig- 
nificant for cold-atom systems. 




FIG. 2. The scaled condensation energy AE/Q in units of 
akp as a function of the filling factor S — N/Q with $1 = 50 
(red squares) and = 200 (blue circles) calculated from exact 
solutions. The scaled interaction is fixed with VqQ = 0.5. The 
solid black curve is the results from the MFT. 

At this stage, we generalize the drifting effect of the 
Zeeman field to the SO-coupled spinor Bose-Einstein con- 
densates, considering that the SO coupling has been real- 
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ized in the ultracold bosonic gas lately The Hamil- 
tonian of the interacting Boson gas with SO and Zee- 
man couplings is taken to be the same as Eq. ((T]) except 
that Ci(;^(^)'s are now annihilation operators of bosons and 
Vo ^ indicating a repulsive interaction between bosons. 
The single-particle spectrum consists of two branches as 
given by Eq. ([6]) with Ck = k-^/2m. According to the dis- 
persion £k,±, one can readily find that the ground state 
of the non-interacting (Vb = 0) ideal boson gas corre- 
sponds to the condensation of all bosons in the lowest- 
lying single-particle state. This energy minimum locates 
at kg — moB in the lower branch (s = —1) with B de- 
notes the direction of the Zeeman field. Therefore, the 
role of the Zeeman field is to lock the direction of the mo- 
mentum of the condensate. Considering the condensation 
of bosons, we can truncate the interaction Hamiltonian 
and obtain the exactly solvable pairing Hamiltonian for 
2D bosons |2d-l23l, which reads 



H= J2 ^k,s<3ak.s + "^^o^4(koMk'(ko), (17) 

k,s=± k,k' 



with s-wave inter-branch bosonic pairing 



%+ko, + °'-k+ko, 



which is distinct from the chiral p-wave intra-branch pair- 
ing of the fermionic case, due to different commutation 
relations. This pairing Hamiltonian takes into account 
the effect of the interaction on the depletion of BEC own- 
ing to the inter-branch pair scattering process of bosons 
from the lowest-lying ko state to the excited states. Note 
that the Hamiltonian is still exact solvable even after in- 
cluding the pair scattering. According to Richardson's 
ansatz [2lj, the eigenstates take the form as 



the fermion (boson) case. This difference originates from 
the important sign change in the commutation relations, 
namely [Ak(ko), (ko)] = S^M'i^ + Es=± "sk-Hko,s] for 
the boson case, in contrast to Eq. pT|) for the fermion 
case. As for the similar current-carrying ground state, 
the underlying physics for both cases is also quite dif- 
ferent. For the fermion case, the fermiology governs the 
ground state (at least in the weak-coupling limit) and the 
nesting of the shifted and/or deformed Fermi surface by 
the Zeeman field favors the Cooper pairing with nonzero 
center-of-mass momentum; while for the boson case, the 
physics is dominated by the tendency of occupation of all 
bosons on the lowest-lying state whose position is deter- 
mined by the Zeeman field. 

To conclude, we have studied the BCS-type pairing 
model of Zeeman-coupled Dirac-type fermions in three 
dimensions and found that the model is exactly solvable. 
In particular, we have revealed rigorously an unconven- 
tional FFLO ground state with the translational momen- 
tum of the Cooper pairs oriented in the direction of the 
Zeeman field. An analogue ground state has also been 
studied for the SO- and Zecman-coupled Bose-Einstein 
condensate. 
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KWk,.})-n^'^n(^ 



t 

k+ko 



k.s 



where 



k ^k-|-ko,+ 



4(ko) 



-k-hko 



(18) 



(19) 



and TOk,s denotes the number of unpaired bosons. Here 
Eij^s obey the following Richardson's equation, 



l + '^c 



1 + mu,+ + "^-k,- 

7b ~7b r 

k "k+ko,+ ^ ^-k-Hko,- 



2Vo 



~^ En — E[y 



= 0. 



(20) 

Despite the seemingly analogy of the pairing Hamil- 
tonian (|8]) and (|17p as well as the resulting Richard- 
son's equations ([M]) and (|20p . some significant differ- 
ences should be noted between the fermion and boson 
systems. The summation over k in the second term of the 
Richardson's equation excludes (includes) the unpaired 
states and the sign preceding the third term is — {+) for 
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